The Markovian metamorphosis of a simple turbulent cascade model 
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Markovian properties of a discrete random multiplicative cascade model of log-normal type 
are discussed. After taking small-scale resummation and breaking of the ultrametric hierarchy 
into account, qualitative agreement with Kramers-Moyal coefficients, recently deduced from a fully 
developed turbulent flow, is achieved. 
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Phenomenological modelling of the energy cascade in fully developed turbulence has a long tradition [jU^}. As 
representatives of the multifractal approach random multiplicative branching processes [^-^| mimic the redistribution 
of energy flux from the large, integral length scale L down to the small, dissipative scale rj and focus on the scaling 
aspect of the surrogate energy dissipation field, extracted from measured velocity time series. In a particular simple 
model version, which is only one-dimensional and binary discrete, a domain of length rj = L/2 J is split into two 
subdomains of equal length r J+ i and the energy flux density e(rj) of the parent domain is non-uniformly redistributed 
by assigning a left/right multiplicative weight qc/n to the left/right subdomain: £c{ r j+i) = qc£(rj), £n{ r j+i) = 
qtze(rj). The multiplicative weights are drawn from a scale- independent symmetric probabilistic splitting function 
p{<lc,qn) = p(qiz,qc), are (qc) = (liz) = 1 on average and are completely uncorrelated to multiplicative weights 
from all other branchings differing in scale and position. At an intermediate length scale r]<rj<L, corresponding to 
j cascade steps, the local bare field density e{rj) = q\q<i ■ ■ ■ qj is a product of j independent multiplicative weights, 
where e(ro) = 1 has been chosen for simplicity. Upon taking the logarithm, the product turns into a summation over 
independent and identically distributed random variables: 

In e(rj ) = In q\ + In q^ + . . . + In qj . (1) 

Introducing yj — y(lj) = lne(rj) — (lne(rj)) as new field variable and lj = \n(L/rj) = jln2 as a logarithmic scale, it 
is straightforward to derive 



Ay Vj+i - Vj 1 , , / (In 2 q) - (lng) 2 

a7 = l^J- = LT2 (ln * +1 - (ln9>) = V 2*2 (2) 

from (Q) for parent/daughter field variables. The last step, leading to the stationary Gaussian- white noise "random 
force" £j with normalisation (£j£j>) — f^^jj' @, only holds for a splitting function p(qc,qTi) — p(qc)p(qTz) of 
log-normal type, where 

K?) = exp I - — ( li i ,/ - — ) j • ( >! ) 




The Langevin equation (g) represents a discrete Markov process, evolving from large to small scales with zero drift 
term = and constant diffusion term — ((In 2 q) — (lng) 2 )/21n2. 

In several ways this thinking seemingly contradicts the results deduced from a large-Reynolds number helium jet 
experiment (7]]^] : from the coarse-grained one-dimensional surrogate energy dissipation field 

1 r x +% 

e(x,r) = - e(x ,r])dx , (4) 

r .L_r 



entering into the transformed variable 



the Kramers-Moyal coefficients 



y(x, I) = lne(x,r) - (lne(x, r)> , (5) 



D^(y(l)) = lim r~ry / ^(2 + M) -y(l)] n p(y(l + M)\y(l))dy(l + Al) (6) 
have been determined to yield 

D^(y)= 7 y (7 « 0.21) 

D {2) {y) = D (Dw0.03) (7) 

D {n - 3) (y) w . 

This outcome suggests that the energy cascade in fully developed turbulence can be described by a scale-continuous 
Markovian Ornstein-Uhlenbeck process, which differs from the discrete random multiplicative branching picture in 
two ways: scale-continuous evolution with a linear drift and constant diffusion term as opposed to a scale-discrete 
evolution with a zero drift and constant diffusion term. - However, the comparison between these two, apparently 
contradicting pictures is not well taken and is not as straightforward as anticipated. So far it has been like comparing 
a caterpillar with a butterfly. Now, we will initiate the metamorphosis in two steps. 
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Step one has to do with the distinction between a bare and a dressed field [Ql . Since fully developed turbulence is a 
three-dimensional process, the redistribution of energy flux from larger to smaller scales should be conserved in three 
dimensions, as long as the dissipative scale rj is not reached. Of course this does not hold once the process is looked 
at in only one dimension, which is done in one-point time-series measurements of one component of the velocity field 
and from which the one-dimensional surrogate energy dissipation field is extracted. For the simple binary discrete 
random multiplicative branching process this implies that the splitting function p{qc, qn) ~ p{qc)p{qn) more or less 
factorises, where p(q) should be a positively skewed distribution limited to the support 0<g<(7 max with g max ~ 2 
p0| . In this respect, the log-normal distribution (||) with the realistic parameter a = 0.42, reproducing observed 
lowest-order scaling exponents (e n (rj))^(L/rj) T ^ and observed multiplier distributions, represents a fair candidate. 
With a factoriscd splitting function, where (qc+qiz) = 2 only holds on average, we have to distinguish between the 
bare field ([!]), which is evolved from the large scale L down to the intermediate scale rj, and the dressed field 

g(r J )=e(r i )(l + A(r i )) ) (8) 

which has been evolved from L all the way down to the dissipative scale rj and then again resummed up to the 
intermediate scale rj. The two fields differ by a small-scale resummation factor (1 + A(rj)), which is equal to one 
only on average. As the experimental analysis (Q) corresponds to the dressed field, we also need to employ the dressed 
field for the random multiplicative branching process, in order to make a fairer comparison between model and data 
results. 

For a truly fair comparison between model and data results we have to call for an additional, second step: since 
the binary discrete random multiplicative branching process is organised hierarchically in one-dimensional space, the 
underlying ultrametric does not allow for spatially homogeneous observables right away. This has been noted only 
recently ]lO|-[T^| and a simple scheme has been suggested, breaking the ultrametric hierarchy and restoring spatial 
homogeneity. It builds a long chain of independent cascade field realisations, each of length L, randomly places the 
observational interval of length r\<r<L within this chain and samples over these random placings. 

These two steps have been decisive for the correct interpretation of the observed multiplier phenomenology 

p3| , pT| : the small-scale resummation of step one explains the scale-independent multiplier distributions as fixed-point 
distributions and step two is in charge for producing the correct correlations between multipliers. - Since the Kramers- 
Moyal coefficients (0) can be understood as moments of logarithmic multipliers, defined for an infinitesimal scale step, 
we might already begin to speculate here that steps one and two might also be in charge for turning the caterpillar 
(g) into the butterfly (g). 

We will now test this speculation by numerical simulation of the binary discrete random multiplicative branching 
process with the factorised splitting function (^) of log-normal type (a = 0.42). A chain of Nl = 10 6 independent 
cascade realisations is constructed, where each realisation has been obtained after J = 10 binary cascade steps; 
consequently the length of the total chain amounts to 10 6 L = 1.28T0 9 7y. 

At first we test the Markov property in general. The conditional probability distributions p(jj(j2)\y{h)) with 
centred intervals h > h is sampled over e "-Nl random placings within the long chain of cascade realisations. It 
is found that these conditional probability distributions fulfill the Chapman-Kolmogorov equation p{y{h)\y(h)) = 
J P^{h)\y{h))p(jj{h)\y{li))dy{h) close to perfectly in the scale range rj<^l<L. This is a necessary and almost sufficient 
validation for this branching process to appear Markovian H . 

Without any loss of generality the Kramers-Moyal coefficients (^) are only calculated at binary scales rj , but again 
sampled over randomly chosen x- values within the long cascade chain. Convergence is tested by letting the positive 
integer number m— >1 in the centred daughter interval of length rj—Arj = rj—2mrj. Good convergence is achieved for 
0<j<5, i.e. the upper part of the cascade inertial range, whereas for 6<j< J = 10 convergence has been found to be 
unsatisfactory since rj-^rj is no longer fulfilled. Consequently, in the following we will only show results for the former 
scale range. - The first Kramers-Moyal coefficient (y, I) is illustrated in Fig. la at binary scales j = 1,3,5. It is not 
constant zero anymore; now it is linear in y(lj). Fitting the parametrisation D^'(y, I) = 7o(0 + 7(01/ yields jo(l) = 
within simulation error bars and a positive, slightly scale-dependent drift coefficient with values listed in Tab. 
la. The latter decreases from a value 0.19 at L = 2 9 i] to a value 0.08 at r$ = 2 5 7] and agrees with the experimentally 
deduced value ^ within a factor of 1.1-2.5. The second Kramers-Moyal coefficient D^ 2 \y, I), depicted in Fig. lb, turns 
out to be almost constant and almost scale-independent; fitted values of the parametrisation (y, /) = D(l)+di(l)y 
can be found in Tab. la. Compared with D = 0.127 of the caterpillar thinking (|^) it is reduced by about a factor 
of 3, but is still about a factor of 1.4 above the experimentally given result (g). Also in qualitative agreement with 
the experimental results, higher-order Kramers-Moyal coefficients are close to zero: 

£> (3) (lM) and£>( 4 %,0 are of the 

order 10~ 3 and 10 -4 , respectively. Here we might evoke Pawulas theorem M to conclude that, after taking steps one 
and two into account, the binary discrete random multiplicative branching process of log-normal type appears as a 
scale-continuous Markovian Ornstein-Uhlenbeck process with Kramers-Moyal coefficients 
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£ (1) (y) = 7(0y 

D [2 \y) = D 
Of" y »(j)«0. 



(0.2 > 7 (0 > 0.1) 
(D re 0.04) 



(9) 



This result is in nice qualitative agreement with the experimental observation 

The result (|^) has been obtained by taking both steps, one and two, of the metamorphosis into account. Step one 
with its small-scale resummation is definitely responsible for the transition from the scale-discrete evolution of the 
bare field to the scale-continuous Markov description for the dressed field. What are then the implications of step two, 
i.e. the breaking of the ultrametric cascade hierarchy to restore spatial homogeneity? In order to clarify this point, we 
restrict the sampling of the Kramers-Moyal coefficients only to the hierarchical positions x m = (m + 0.5)rj with integer 
Q<m<2 3 ■ Nl within the long chain of cascade configurations. For these positions, the integration interval of length rj, 
entering into (Q), perfectly matches an intermediate interval of the bare cascade evolution. Results for the first and 
second Kramers-Moyal coefficients are listed in Tab. lb. The first coefficient, which we now denote with a tilde, is again 
found to be of the form D^(y,l) — j(l)y- Note however, that the drift coefficient 7 is negative and that its modulus 
is about a factor of 3-7 less when compared with 7. This demonstrates that small-scale resummation alone already 
introduces a weak linear drift term, but breaking of the ultrametric hierarchy is very necessary to change its sign and 
to bring it to the correct order of magnitude. Also the second Kramers-Moyal coefficient D^ 2 \y 7 l) = D(l) + di(l)y 
is affected by leaving out step two, but only weakly: again it shows a small scale-dependence, is almost constant for 
a given scale and, once compared with D(l), is reduced by about a factor of 1.2-1.7. Note also, that it is mainly 
small-scale resummation, which drives the diffusion coefficient away from the caterpillar thinking (||) with D = 0.127. 

We conclude: small-scale resummation and breaking of the ultrametric hierarchy initiate the Markovian meta- 
morphosis of a discrete random multiplicative branching process as they turn the caterpillar, a discrete, Gaussian 
white-noise evolution in scale, into a butterfly, that is an effective scale-continuous Ornstein-Uhlenbeck description, 
the latter being in qualitative agreement with the experimental observations. At least on a qualitative level, there is 
no conflict between the experimentally observed energy cascade in fully developed turbulence and random multiplica- 
tive branching processes; this statement is further supported by recent work on multiplier distributions ]Tc|Jl2] , ^5[ . 
- Several points have to be considered in order to possibly achieve an even better quantitative agreement between 
such model results and experimental observations: sensitivity on the generator of random multiplicative branching 
processes, i.e. on scale-discrete or scale-continuous implementations and on the choice of splitting function, as well 
as on finite-size, large-scale and dissipative effects; work in these directions is in progress |l6|. Also data results for 
larger Reynolds numbers would be highly appreciated. 
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di 
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di 





0.154 


0.042 


-0.004 


-0.178 


0.036 


0.004 


1 


0.193 


0.042 


-0.004 


-0.077 


0.034 


0.002 


2 


0.163 


0.040 


-0.003 


-0.052 


0.031 


0.001 


3 


0.133 


0.040 


-0.002 


-0.034 


0.029 


0.001 


4 


0.105 


0.040 


-0.002 


-0.021 


0.027 


0.001 


5 


0.084 


0.043 


-0.002 


-0.012 


0.025 


0.001 



TABLE I. First and second Kramers-Moyal coefficients, D {1) (y,lj) = j(lj)y and D^'iyjj) = D(lj) + ^(l^y, after 
small-scale resummation and (a) with / (b) without breaking of the ultrametric hierarchy. 
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FIG. 1. First (a) and second (b) Kramers- Moyal coefficients, obtained after small-scale resummation and breaking of the 
ultrametric hierarchy. The Z)' 2 ' -values for j = 3 and 5 have been artificially shifted by 0.02 and 0.04, respectively. 
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